A Quantum Optical Spring 
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We study the dynamics of the quantum optical spring, i.e., a spring whose spring constant under- 
goes discreet jumps depending on the quantum state of another system. We show the existence of 
revivals and fractional revivals in the quantum dynamics reminiscent of similar dynamical features 
in cavity QED. We recover in the semi classical limit the results for an oscillator whose frequency 
undergoes a sudden change. The quantum optical spring is conceivable for example by a micromirror 
under the influence of radiation pressure by a field which is strictly quantum. Our work suggests 
that driven systems would in general exhibit a very different dynamics if the drive is replaced by a 
quantum source. 

PACS numbers: 42.50.-p, 42.50.Pq 



I. INTRODUCTION 

The quantum dynamics of an oscillator with time de- 
pendent modulation of its frequency has been investi- 
gated very extensively [J, H, 0, H[ . This is because of its 
many applications in different contexts. One very inter- 
esting application was in the context of center of mass 
motion of an ion in Paul trap. More recently the study 
of modulated oscillator acquired new importance in con- 
nection with BEC's where trap potential can easily be 
manipulated 0, Q. Further nanomechanical oscillators 
provided new impetus for the study of quantum dynamics 
of frequency modulated system Q • One of the important 
outcome of such studies was that mechanical and vibra- 
tional degrees of freedom can be prepared in nonclassical 
states In particular, one found that sudden changes 

in the frequency of the oscillator yield significant amount 
of squeezing. 

In all the above studies the oscillator was treated as a 
quantum system whereas the source of modulation was 
considered as something prescribed. Clearly a complete 
quantum dynamical theory should consider the source 
also to be quantized system. This is important as the 
mechanical oscillator can have significant back action on 
the source of modulation. The system is then a quan- 
tum optical spring , i.e., a spring whose spring constant 
depends on the quantum state of another system. The 
modulation by a quantized source produces new features 
in the quantum dynamics as the oscillator samples the 
discrete structure of the quantized source of modulation. 
In particular we show the appearance of characteristic 
quantum revivals and fractional revivals 0, [l^, [H| ■ The 
dynamics we present should be realized by the effects of 
radiation pressure on micromirror 

The organization of the paper is as follows : In section 
II, we describe the model and derive analytical result 
for the reduced density matrix of the oscillator and the 
source of modulation. The quantum dynamics of the 
oscillator is discussed in section III. In section IV, we 
analyze the squeezing properties. The back action effect 
is described in section V. 



II. MODEL FOR THE QUANTUM OPTICAL 
SPRING 



Consider a harmonic oscillator with mass m and fre- 
quency u). Its Hamiltonian is given by 



Hn = 



1 9 9 

p mu> x 
2^ + 2 



(1) 



where x and p are the position and momentum co- 
ordinates obeying the canonical commutation relation 
[x,p] — ih . Let us assume that the frequency of the 
oscillator is modulated by a quantized source such that 
the Hamiltonian in Eq. |T]) goes to 



H 



p 2 mw 2 (l + iia^a)x 2 
2^ + 2 



(2) 



Since the interaction part commutes with the unper- 
turbed Hamiltonian of the quantized source we do not 
write this part of the energy. Further for brevity we 
would denote by QSM as the quantized source of mod- 
ulation. The operators a and for the QSM obey the 
Boson commutation relations [a, a^] = 1. The Hamilto- 
nian @ describes the quantum optical spring. Let the 
eigenstates and eigenvalues of H be denoted by 



4> n {x) = N n H n (ax)exp (~-a 2 x 2 ), 



a 

Vw2"n! ■ 



/ a n' 



(3) 



where H n is the nth Hermite polynomial. 
The eigenstates of H are given by 
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a T a|p) =p\p), 

uj p = u>y/{l + lip) . (4) 
and where the wave function ip? is given by 



fa 1 

i>n = N n H n {a p x){J — )exp(--a 2 :r 2 ), 
V a & 

Note that ip? is the eigenstate of harmonic oscillator with 
frequency to replaced by ui p . For a fix p, these form a 
complete set. 

The time evolution of the initial state can be obtained 
from the knowledge of the states given by Eq. Q. Let 
us consider the initial state of the coupled system given 
by 



ip(t = 0) = Y, C p^{x)\p) , 



(6) 



where C pn are the expansion coefficients. The state at 
time t would then be 



\m) 



exp(— — )xP(t = 0) 



~ C pn \p) exp[— (|- + —?—)}cp n (x) 



^C p „|p)exp[— ( 



h 2m 



? 2 

mujpX , 



—itF p 
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(7) 



This is the state of the combined system of oscillator 
and the QSM. It is a nonfactorized state and thus the 
spring gets entangled to the quantum system controlling 
the spring constant. The reduced state of the oscillator 
and QSM can be obtained by projecting out the degrees 
of freedom of the other system. It is seen that the density 
matrix for the oscillator is 



p = £ m(^\c pn c; r , 



x exp[- 



-it(Ef - E p .) 

^}(^\<p n )(4> m m , (8) 



and that for QSM is 



PQSM = ^ d PiP2\Pl)(Pl\ > 

Pl,P2 



(9) 



^Plp2 ~~ Q>inCp 2m ('0^ 2 |' i / 7 f 1 / 



X^^nX^mlVf >exp 



-it{E^-Ef ) 



(10) 



Note specially the rather involved form of the reduced 
density matrix for the QSM. This is because the scalar 
product of the wave functions ipf for different p is 
nonzero- 



(Vf ¥= for p 1 ^p 2 

= 5i j if pi = p 2 



(11) 



Using Eq. ([8]) we can study the details of the quantum 
dynamics of the oscillator coupled to QSM. We would like 
to mention that the quantum characteristics of a nanome- 
chanical oscillator coupled to a quantized photon source 
has been studied by several authors [3, [lU, [H, 03] who 
have shown how nonclassical states like CAT states of 
such a system can be generated. The Hamiltonian used 
by previous authors is given by a) ax which is linear in 
the oscillator variable and is different from a^ax 2 which 
is produced by modulation. However Arcizet et al. fl2j 
have used precisely this Hamiltonian as they argue that 
the radiation pressure leads to modulation of the spring 
constant. Another possibility to realize quantum optical 
spring would be the center of mass motion in a detuned 
cavity. 



III. QUANTUM DYNAMICS OF THE 
QUANTUM OPTICAL SPRING : REVIVALS 
AND FRACTIONAL REVIVALS 

Let us consider first a simple case when the QSM is 
prepared in a coherent state with mean amplitude a and 
the oscillator is prepared in its ground state. Thus the 



coefficients C pn are given by 



a 



pa 



I r 2 

a p exp( L i ) 



(12) 



Using Eq. (fT2"j) in Eq. (JSJ) and on simplification we get 



Po = 



E 



a| 2 Pexp(-|a| 2 ) 



(P)l 



■expl 



x(^o}(0o|^> 



(13) 
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And in particular the probability of finding the oscillator 
in the initial state is given by 



P (t) = (Mpo\M 



|a| 2 Pexp(-|a| 2 ) 



El ^ 



(P)! 

2 |a|^exp(-|a| 2 ) 



(14) 



where 



A> = ^2exp[-itui p l}\ (ipf\ 



(15) 



This is one of our key results. The quantum state of 
the source appears through the weight factor in Eq. (fT4|) 
and through the discrete sum in Eq. (|14p . Such sums 
have been encountered before in connection with the cav- 
ity QED problems d G3 and have been experimentally 
studied [H, 0j[ . The sum in Eq. (fT5|) can be evaluated 
in closed form since 



(Vf |<£o) = / ip?(x)<l>o(x)dx 



x / H ; (apir)exp[(--(a 2 + a 2 )x 2 )]cfo 



NtN (-&=) 

x y Hi(f3 p y)exp[-y 2 }dy 

JOLpOL [my. 



f3 p VW(^-1) T 



where 



(r? p )* L 2™(ro)! 

% = (1 + MP) , 
2 2 (%)^ 



(16) 



(17) 



On using Eq. (jl6|) in Eq. (fTSj) we find closed form ex- 
pression for A p : 



\P P \ 2 



^(l-(^-l) a exp(-2iw p t))a 



(18) 
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FIG. 1: Figure shows the variation of Po as a function of 
r(ujt — 2nr) with parameter /i = 0.1. The average number 
| ct | 2 of excitation in the QSM is 4. 
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FIG. 2: Figure shows the variation of Po as a function of 
r(ujt = 27tt) with parameter /j, = 0.3. The average number 
|q| 2 of excitation in the QSM is 25. 



and Pq(£) can be written as 



Po(t) = (\A P \ 2 



(19) 



and where () is the average over the distribution of p. 
For the classical source of modulation the corresponding 
result would be \A p \ 2 with p replaced by the strength of 
modulation \a\ 2 



Pd{t) 

(3 a 2 



\(3 a \ 



y/ruR - 2(/? 2 - l) 2 cos(2^i) + (/? 2 - l) 4 ) 
(l + /x|a| 2 ) , 



2(»?a) 



uo a = toy/ (1 + n\a\ 2 ) 



(20) 



Thus P c i oscillates at the frequency 2 uj a . However the 
probability Po(t) for the QSM behaves quite differently 
as shown in the Fig. 1 and exhibits collapse and revival 
of the classical periodic motion. In a different parameter 
regime the fractional revivals in Po(t) are seen in Fig. 2. 
Clearly the dynamics of the quantum optical spring is 
very sensitive to the quantized source of modulation. 
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IV. SQUEEZING PROPERTIES OF QUANTUM 
OPTICAL SPRING 



We next investigate if the oscillator system acquires 
important nonclassical character as a result of QSM. As 
mentioned earlier even a classical change in the frequency 
can lead to squeezing in oscillator and therefore an inter- 
esting question would be how this squeezing character is 
modified by the QSM. In order to calculate squeezing it 
is more convenient to work with Heisenberg equations of 
motion. For the Hamiltonian in Eq. (f5]) we can prove 
that 
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FIG. 3: Figure shows the variation of V x (t) as a function of 
r(ut = 2ttt) with parameter fi = 0.1. The average number 



of excitation in the QSM is 4. 



d(a^a) 
dt 

dx 

~dl ~ 



. 



P_ 

m 



= — m(cj) 2 (l + [ia)a)x , 



(21) 



(22) 



(23) 



Thus a^a is a constant of motion and we can write the 
solution as 
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FIG. 4: Figure shows the variation of V x (t) as a function of 
T(tut = 2ht) with parameter fi = 0.3. The average number 
\a\ of excitation in the QSM is 25. 



x(t) = ie(0) cos(Of) + sin(fit) , 



p(t) = p(0) cos(Qi) - mCl x(0) sin(fii) 
where is given by 



Cl = cj\/l + ua+a 



and where a^a is the operator at time t = 0. 
We define the squeezing variables 



V x (t) = 



v P (t) 



{x 2 {t)) (x(t)) 2 

<z 2 (0)} 
(P 2 (t))-(p(t)f 



(P 2 (0)) 

For the ground state of the oscillator 



(x(Q)) = (p(0)> = 0,(xp+px) = 0, 
(* 2 (0)> 



2mw 

Using Eqs. ((Ml) and (TJ7D we get 



(p 2 (0)) = -s- 



(24) 



(25) 



V x (t) = 1- 



' exp[— |a| 2 ] sin [wt-^/l + /m ] /in 



fin 
(28) 



v P (t) 



y-v |a| 2 "exp[ 



|a| 2 ] sin [wt-^/l + un ](un) 



(29) 

Clearly V x (t) is always less than one and hence the x 
quadrature is squeezed. The corresponding result with a 
classical source of modulation would be 



V x {t) = l-sinWl + MM 2 K 



1 + H\a\ 



(30) 



The minimum value of V x in classical case is when 



(26) sin [uity/1 + fi\a\ 2 ] = 1 and hence 

x ' mm l l + u|a| 2 



)■ 



(31) 



For the QSM the squeezing parameter (|2"8|) exhibits typi- 
cal collapse and revival of the classical periodical oscilla- 
(27) ti° n Pigs. 3, 4]. This is again due to the discrete nature 
of the quantum state of the source of modulation, i.e., 
control of the spring constant by a quantum source. 
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V. BACK REACTION OF THE OSCILLATOR 
ON QSM 

In this section we evaluate the back action of the oscil- 
lator on the QSM. The Hamiltonian in Eq. @ depends 
only on a^a and therefore only the off diagonal elements 
of the density matrix pqsm are affected by the coupling 
of the oscillator to QSM. The reduced state of the QSM 
is rather involved. In order to appreciate the back ac- 
tion we examine the mean displacement (a) of the QSM 
which is conditional on the measurement of oscillator in 
the state \4>q). The mean value is given by 



d(t) ee (a) 

= £<J|a|n>(n|p|0 

1 1. 1 

= ^ VnPn,n-l , (32) 
n 

where QSM's conditional density matrix is found to be 

Pn,l = ^ c„ciX nt i , (33) 

with X n> i given by 

X n>i = ((p \exp[ih(n)]exp[-ih(l)]\4>o) . (34) 

In the above equation we have introduced the Hamilto- 
nian h(n) defined by 



h(n) 



h 2m 



mw 2 (l + {j,n)x* 



) 



(35) 



The X n ys can be calculated by using operator disentan- 
gling theorems for SU(1,1) group . We give details in the 
Appendix(A). We find the following result [Eq. (|A13j) ] 



X n t = 



((r 3 n)*r 3i )4 
- (r+„)*r 



(36) 



+/ 



The displacement of the QSM with Ps defined by Eqs. 
(|A9|l and (fAlO)) can be obtained from Eqs. ([3"2]) , (jS3"]) . 

In Fig. 5 we display the time evolution of the imag- 
inary part of (a). Note that the imaginary part of (a) 
builds up starting from zero and exhibits aperiodic be- 
havior. The quantized nature of the optical spring is well 
reflected in the back reaction [20] . Evidently the back re- 
action occurs over a much larger time scale. 

In conclusion we introduced the idea of a quantum op- 
tical spring , i.e., a spring whose spring constant is con- 
trolled by another quantum system. The dynamics of the 



oscillator system exhibits the phenomena of collapse and 
revivals including fractional revivals. Further the spring 
gets entangled to the quantum system which controls the 
spring constant. 
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FIG. 5: Figure shows the variation of Imaginary (a) as a func- 
tion of r(u)t = 167tt) with parameter [s, — 0.3. The average 
number \a\ 2 of excitation in the QSM is 25. 



APPENDIX A: CALCULATION OF THE 
MATRIX ELEMENT (f34|) 

In the appendix we give the details of our calculation 
for X n i- It is convenient to use x and p in terms of 
angulation and creation operators b and w 



2muj 



(Al) 



expHM«o) = ex P [-f (|1 + ^ 2(1 + ^ 2 )]|<M 



(A2) 



Substituting the values for x and p from Eq. (jAIj) in the 
above equation and on simplification we get 



exp[-i/i(n)]|^ ) = exp[y +n K + +j- n K- +l3nK 3 }\4) Q } 



(A3) 



where 



7±n 



73n 



-ipumt 



-iuit(p,n + 2) 



(A4) 



(A5) 
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In (|A3[) K±, K 3 are the generators of SU(1,1) algebra given by 

K 3 = ±@b + ffi) , (A6) 

K+ = {k^ = \{tfY . (A7) 
The Eq. (|A3|) can be simplified by using the disentangling theorem 21] for the SU(1,1) group: 

exp( 7+n X + + 7_ n #_ + i 3n K 3 ) = exp(r +n if + ) exp[(ln T 3n )K 3 ] exp(r_„if_) , (A8) 

where 

F3 " = (cosh/3„-^sinh/3„)2 ' (A9) 

r±„ ee 2 7±ra sinh/3» ^ 

(2/3„ cosh/3„ - 73„ sinh/3„) 

^n 2 = ^l3n 2 ~ 1+nl-n • (All) 

We use the disentangling theorem in the form (|A8[) as then acting on |</>o) yields zero and K 3 terms can be 
simplified as |0o) is an eigenstate of K 3 . Thus we reduce Eq. (IA3|) to 

exp(j +n k + + j- n K- + 73»j^3)|0o) = exp(r + „X + ) cxp[(ln r 3n )A" 3 ] exp(r_ n if_)|<?io) 

= exp[^^)]exp[r +n ^ + ]|0 o ) 

= (T 3n y exp(r+„K + )|0 o ) (A12) 
Hence X n j as defined by ([34]) can be calculated as follows 



X„j = (()>o\exp(ih(n))exp(-ih(l))\<f> ) 

= ((r3„)*r 3i )3(0o|exp((r +n )*^_)exp(r +; if + )|0 o ) 

- Ill 3 „j 1 30 4 ( p! )2 2 2 P 

- « r »»>' r «>* (A13) 



v/i - (r+„)*r +i 

This is the result we use in Section V to calculate the back-action 
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